Preprint typeset in .IHEP style - HYPER VERSION 



DF/IST-17.2003 
hep-th/0311260 



Scalar-gravitational perturbations and quasinormal modes in the 

^ : five dimensional Schwarzschild black hole 

o 
o 

(N 
o 

. 

Q , Vitor Cardoso^ Jose P. S. Lemos^^ and Shijun Yoshida^^ 
OO ■ ^ 

'Centra Multidisciplinar de Astrofisica - CENTRA, Departamento de Fisica, Instituto Superior Tecnico, 
Av. Rovisco Pais 1, 1096 Lisbon, Portugal 

(N 
> 
O 

Abstract: We calculate the quasinormal modes (QNMs) for gravitational perturbations of the 
^f) _ Schwarzschild black hole in the five dimensional (5D) spacetime with a continued fraction method. For 
O ■ all the types of perturbations (scalar-gravitational, vector-gravitational, and tensor-gravitational pertur- 
. bations), the QNMs associated with / = 2, / = 3, and / = 4 are calculated. Our numerical results 
I ' are summarized as follows: (i) The three types of gravitational perturbations associated with the same 
angular quantum number / have a different set of the quasinormal (QN) frequencies; (ii) There is no 
> purely imaginary frequency mode; (iii) The three types of gravitational perturbations have the same 
J> ' asymptotic behavior of the QNMs in the limit of the large imaginary frequencies, which are given by 
'Cl ; ^'^H^ ^ logs 27ri{n + 1/2) as n ^ oo, where w, Th, and n are the oscillation frequency, the Hawking 
■ temperature of the black hole, and the mode number, respectively. 



vcardosoOf isica. ist.utl.pt, LemosQkelvin.. ist.utl.pt, yoshidaOf isica. ist . utl .pt 



Keywords: Quasinormal Modes, Gravitational Radiation, Higher Dimensions 



Contents 

|l|. Introduction 

^. Numerical Method 

^. Numerical Results 

^. Conclusion 



1. Introduction 



The quasinormal (QN) ringing of a Schwarzschild spacetime was first observed in through numerical 
calculations of the gravitational wave scattering by the black hole. Since then, the quasinormal modes 
(QNMs) of black holes have been extensively studied. The classical motivation behind the exploration of 
the QNMs of black holes is twofold: One is to answer the question of whether the spacetime is stable, and 
the other to know what kind of oscillations will be excited in the spacetime as some perturbations are given. 
The latter is quite important from the observational point of view because we could determine fundamental 
parameters of a black hole, such as the mass or the angular momentum, through the information of the 
QNMs. A number of studies on QNMs of several different spacetimes containing black holes have been done, 
(for a review see, e.g., [@, Recently QNMs have acquired a different status, since it was conjectured 
that they may be connected to black hole area quantization and quantum gravity @, 

Most studies on QNMs of black holes were restricted to the four dimensional (4D) case, compatible 
with astrophysical scenarios. However, motivated by the TeV-scale gravity proposals for instance, 
higher-dimensional (higher-D) theories of gravity have recently attracted much attention. Up to recently, 
there were no master equations for examining the QNMs for gravitational perturbations of higher-D black 
holes. Due to the absence of these master equations, only the QNMs of the simplest test fields, namely 



massless scalar fields, around the higher-D black holes were calculated [^, ^, 11 1. 

The situation has changed recently. Kodama and Ishibashi jl^ have derived master equations for grav- 
itational perturbations in a higher-D Schwarzschild black hole spacetime. Putting D = 2 + n, it was shown 
that for the case n > 3, the gravitational perturbations can be divided into three classes, namely, scalar- 
gravitational, vector-gravitational, and tensor-gravitational perturbations, according to their tensorial be- 
havior on the n-sphere. The scalar-gravitational and vector-gravitational perturbations correspond to the 
polar and axial perturbations in the 4D spacetime. The tensor-gravitational perturbations are a new kind 
appearing in higher-Ds. Kodama and Ishibashi have also shown that all the types of gravitational perturba- 
tions can be reduced to a simple Schrodinger-type wave equation like Regge- Wheeler or Zerilli equations. It 
is important to ask the question of whether there is a special relationship among the scalar-gravitational, 
vector-gravitational, and tensor-gravitational perturbations because in 4D the scalar-gravitational and 
vector-gravitational perturbations have a special relation yielding the same quasinormal (QN) frequencies. 
There seems to be no such relation between the scalar-gravitational and vector-gravitational perturbations 
11. 

Recently, some studies on the QNMs for gravitational perturbations of higher-D Schwarzschild black 
holes have appeared. Konoplya has calculated the fundamental QNMs of Kodama and Ishibashi's master 
equations with a JWKB approximation, and has found that three types of perturbations have different 
fundamental QN frequencies [p^ j. Berti, Cavaglia, and Gualtieri have also done similar calculations but for 
a wide range of angular eigenvalues I of perturbations [0]. Cardoso, Lemos, and Yoshida have calculated 
the QNMs for the vector-gravitational and tensor-gravitational perturbations up to higher-order modes 
1^]. As for the asymptotic behavior in the limit of highly damped modes, it has been shown that the three 
types of gravitational perturbations have the same behavior regardless of the spacetime dimensions or the 
angular quantum numbers of the perturbations |l^, ^. 

The purpose of this paper is to go a step further up and explore the QNMs in the special 5D higher-D 
black hole spacetime. In particular, we want to find out whether or not there is a relationship among 
the QNMs of the scalar-gravitational, vector-gravitational and tensor-gravitational perturbations in this 
5D Schwarzschild black hole. The 5D case can be considered representative of the other higher-D cases. 
A prime purpose of the present study is to extensively calculate the QNMs for the scalar-gravitational 
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perturbations. The scalar-gravitational perturbations are the most important because they can be excited 
more easily than other types of perturbations. We calculate the QN frequencies not only for low-order 
modes but also for relatively higher-order modes, enabling us to discuss the asymptotic behavior of the 
highly damped QNMs. In the present study Nollert's approach is employed numerically to obtain the 
QN frequencies. The paper is organized as follows: In section 2 we present the basic equations employed 
for obtaining QNMs in the 5D Schwarzschild spacetime. Numerical results are given in section 3, and in 
section 4 we conclude. 



2. Numerical Method 



The general perturbation equations for the Schwarzschild black hole in D = n + 2 dimensions have been 
recently derived by Kodama and Ishibashi |jl^. It has been shown that there are three completely decoupled 
classes in the perturbations for the n > 3 case, namely scalar-gravitational, vector-gravitational, and tensor- 
gravitational perturbations, and that all types of perturbations can be cast into a Schrodinger-type wave 
equation, given by 

f^(f^)+i-^'-Vir))^ = 0, (2.1) 



dr \ dr ^ 

where the function / is minus the gu component of the metric tensor in Schwarzschild coordinates, and 
given by 

2M , , 

f = ^-^ = ^-^- (2-2) 

Here, r is the radial coordinate of the n + 2 dimensional Schwarzschild spacetime, and uj is the oscillation 
frequency of the perturbations. M is a parameter related to the mass M of the black hole via the relation 

M = -— ^, (2.3) 

ovrGn+2 

where is the area of a unit n-sphere, given by An = 27r*^"'^^)/^/r[(n + l)/2], and Gn+2 stands for the 
n + 2 dimensional Newton constant. The Hawking temperature of the black hole is Th = (n — l)/47rr/i, 
where rh is the horizon radius, defined by r^~^ = 2M. 

In the master equation (p.l|), depending on the type of perturbations, the effective potential can be 
written as 



/ 



Q{r) 



for scalar 



16//(r)2 

(2/+n)(y+n-2) _ Sr^ ^^^^^^ (2.4) 

(2/+n)(2^+n-2) ^ ^ for tcnsor , 



with 



Q(r) = n^{n + ifx^ + n{n + l){4{2n^ - 3n + 4)m + n(n - 2)(n - 4)(n + l)}x^ (2.5) 

— 12n{(n — 4)m + n(n + l)(n — 2)}mx + 16m^ + 4n(n + 2)m^ , 
1 

H{r) = m + —n{n + l)x , (2-6) 

m = l{l + n-l) -n, (2.7) 

where I means the angular quantum number of perturbations, and I > 2 has been assumed. It is worthwhile 

to note that the scalar-gravitational and vector-gravitational perturbations are, respectively, higher-D 
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counterparts of the polar and axial perturbations in the 4D Schwarzschild black hole, while the tensor- 
gravitational perturbations are brand new in the sense that a counterpart of the tensor-gravitational 
perturbations does not appear in the 4D case. In the case of n = 2, the master equations for the scalar- 
gravitational and vector-gravitational perturbations therefore become the Zerilli equation and the Regge- 
Wheeler equation, respectively |17, Interestingly, the effective potential for the tensor-gravitational 



perturbations is exactly the same as that for the massless scalar field in the higher-D Schwarzschild black 
hole |]. 

The QNMs of the higher-D Schwarzschild black hole are characterized by the boundary conditions of 
incoming waves at the black hole horizon and outgoing waves at spatial infinity, written as 

Hr) I Cr " " " (2.8) 

where the time dependence of perturbations has been assumed to be e*"^*. Here, r* denotes the tortoise 
coordinate, defined by dr^ = f^^dr. In order to obtain numerically the QN frequencies, we employ Nollert's 
method [jl8|, an enhanced version of Leaver's continued fraction method since it is nicely suitable to 
the study of the asymptotic behavior of QNMs in the limit of highly damped modes, the main concern in 
the present investigation. 

In this study, we only consider the QNMs of the Schwarzschild black hole in the five dimensional 
spacetime, namely the n = 3 case. The tortoise coordinate is then reduced to 

n = z-^ + ^ ln{z - zi) - -L ln{z + zi) , (2.9) 
2zi lz\ 

where z = r^^ and z\ = r^^. The perturbation function ^ can be expanded around the horizon as 

<I> = e--^"(z-zi)''(z + zi)''^aJ^^ , (2.10) 

fe=o ^ ^ 

where p = iuj/2zi and oq is taken to be oq = 1. Here, the expansion coefficients for k > 1 are 
determined with the recurrence relation, the order of which is dependent on the functional form of the 



eff'ective potential in equation (2.1). As shown by Cardoso, Lemos and Yoshida the recurrence relation 
becomes a four-term relation for vector-gravitational and tensor-gravitational perturbations. In this paper 
we do not show an explicit expression of the four-term recurrence relation for the vector-gravitational and 
tensor-gravitational perturbations since it can be found in |5|. For the scalar-gravitational perturbations, 
on the other hand, the recurrence relation becomes eight-term relation, given by 

co(0)ai + ci(0)ao = 0, 

co(l)a2 + ci(l)ai + C2(l)ao = 0, 

co(2)a3 + ci(2)a2 + C2(2)ai + 03(2)00 = , 

co(3)a4 -I- ci(3)a3 -|- 02(3)02 + 03(3)01 + 04(3)00 = , 

00(4)05 + 01(4)04 + 02(4)03 + 03(4)02 + 04(4)01 + 05(4)00 = , 

oo(5)o6 + 01(5)05 + 02(5)04 + 03(5)03 + 04(5)02 + 05(5)01 + O6(5)oo = 0, 

oo(A:)oa;+i + ci{k)ak + 02(/c)ofc_i + 03(/c)ofc_2 + Ci{k)ak^'i + C5{k)ak-4 

+CG{k)ak-5 + C7{k)ak-e = , for /c = 6, 7 • • • , (2.11) 



-A- 



where 



coik) = 2(1 + k){6 + mf{l + k + 2p), 

ci{k) = -(6 + m){m^ + (78 + 5m)k^ + Wm{l + 2p)p + (30 + 5m + 216/9 + 2Qmp)k + 12(1 + 5p + 10/3^)} , 
C2{k) = [8(324 + 48m + m^)^^ + ^2^1 + 32^2) _^ i2m(13 - 48p + 112^^) + 36(41 - 128p + l^p^) 

+16{-198 + 540p + 2m2p + 3m(-7 + 30p)}A;]/2 , 
C3(/c) = [-4(1980 + 168m + m'^){-2 + kf - m^{'i + 4pf - 12m(69 + 192p + 224p2) 

-36(129 + 520p + 752/9^) - 4(-2 + k){m'^{3 + ip) + 96m(3 + 7p) + 36(65 + 204p)}]/4 , 
c^(k) = 12[6(25 + m)/c2 + m(25 - Up + 24^^) + 6(127 - 2Up + 96^^) 

+{-660 + 588p + m(-22 + 24p)}/c] , 
C5{k) = -3[4(81 + m)A;2 + ^(29 - 40p + 16/3^) + 6(491 - 644p + 216/9^) 

+4{-483 + 324p + m(-5 + 4p)}k] , 
ce{k) = 18{227 + I6k'^ - 240p + 64p^ + 8(-15 + 8p)k} , 
cj{k) = -9{-9 + 2k + 4pf . 

(2.12) 

It is understood that since the asymptotic form of the perturbations as r^, — > cxd is written in terms of the 
variable z as 

g-ic.r. ^ g-i<^.-i _ + Zi)P, (2.13) 



the expanded perturbation function <I> defined by equation ( |2.10 ) automatically satisfy the QNM boundary 
conditions (|2.8D if the power series converges for < z < zi. After making a Gaussian elimination five 



times [20], we can reduce the eight-term recurrence relation to a three-term relation, which has the form 

aoai + Poao = , 

QfcOfc+i + PkCLk + TfcOfc-i = , = 1, 2, • • • . (2.14) 

Here, we omit the explicit expression for the final three-term recurrence relation because the numerical 
procedure to obtain the three-term recurrence relation is quite simple. Now that we have the three- 
term recurrence relation for determining the expansion coefficients a^, according to Leaver [21, 19|, the 
convergence condition for the expansion ( 2.1C| ), namely the QNM conditions, can be written in terms of 
the continued fraction as 

ao7i "172 0273 , 

pi- P2- pa- 
in order to use Nollert's method, with which relatively higher-order QNMs with large imaginary frequencies 
can be obtained, we have to know the asymptotic behavior of at+i/ak in the limit of A; — > cxd. Following 
Leaver [20[, it is found that 



Qfc+i 



l±2^k-^''^ +{2p-^k-^ + ■■■ , (2.16) 
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Table 1: QN frequencies ujrh of the scalar-gravitational perturbations for / = 2, / = 3 and I = 4. 



n 




/ 


= 


2 




/ 


= 


3 






= 


4 








9477 


+ 


0.2561i 


1 


6056 


+ 


0.3110i 


2 


1924 


+ 


0.3293? 


1 





8512 


+ 


0.8212i 


1 


5109 


+ 


0.9528? 


2 


1149 


+ 


0.9999i 


2 





6727 


+ 


1.5431i 


1 


3325 


+ 


1.6586? 


1 


9649 


+ 


1.7078? 


3 





5121 


+ 


2.4399i 


1 


1107 


+ 


2.4693i 


1 


7593 


+ 


2.4804i 


4 





4169 


+ 


3.4140ii 





9091 


+ 


3.3823i 


1 


5319 


+ 


3.3352i 


5 





3619 


+ 


4.4092? 





7586 


+ 


4.3533? 


1 


3215 


+ 


4.2642? 


6 





3275 


+ 


5.4102? 





6524 


+ 


5.3471? 


1 


1486 


+ 


5.2403i 


7 





3043 


+ 


6.4130i 





5765 


+ 


6.3485? 


1 


0130 


+ 


6.2392ii 


8 





2876 


+ 


7.4164? 





5204 


+ 


7.3524? 





9071 


+ 


7.2474i 


9 





2751 


+ 


8.4197? 





4778 


+ 


8.3569? 





8232 


+ 


8.2584i 


10 





2653 


+ 


9.4229i 





4445 


+ 


9.3614? 





7555 


+ 


9.2698i 


11 





2574 


+ 


10.426i 





4178 


+ 


10.366? 





7000 


+ 


10.28H 


12 





2510 


+ 


11.429? 





3961 


+ 


11.370? 





6537 


+ 


11.290? 



where the sign for the second term in the right-hand side is chosen so as to be 

Re(±2^)<0. (2.17) 

In actual numerical computations, it is convenient to solve the k-th inversion of the continue fraction 
equation ( p.l5| ), given by 

afc-i7fc afc-27fc-i "oTi 

Pk 



Pk-i— (3k-2— Po 

QfcTfc+l QA;+l7fc+2 
Pk+l— Pk+2 — 



(2.18) 



The asymptotic form ( p. 161) plays an important role in Nollert's method when the infinite continued fraction 
in the right-hand side of equation ( 2.1^ ) is evaluated |18]. 



3. Numerical Results 

In order to confirm that our approach to obtain the QNMs for the scalar-gravitational perturbations works 
well, we first calculate the QNMs in the 4D case, namely the QNMs of Zerilli equation. For the Zerilli 
equation, eigenfunction is expanded around the horizon as 



$ = e-^'^^"'(z-zi)*^^i"'z^'^^i"' J]]afc (^^] , (3.1) 



k=0 ^ ^ 

and the corresponding recurrence relation becomes a five-term relation. The numerical result obtained are 
in good agreement with the QN frequencies of the Regge- Wheeler equations, which are exactly the same 
as those of the Zerilli equation, showing that our numerical approach works quite well. 

We calculate the QNMs associated with / = 2, 3, 4 for the scalar-gravitational, vector-gravitational, 
and tensor-gravitational perturbations in the five dimensional Schwarzschild black with Nollert's numerical 
approach [08]. In order to discuss the asymptotic behavior of the QNMs in the limit of large imaginary 
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Table 2: QN frequencies ujrh of the vector-gravitational perturbations for 1 = 2,1^3 and I = 4. 



n 




/ 


= 


2 




/ 


= 


3 




/ 


= 


4 





1 


1340 


+ 


0.3275i 


1 


7254 


+ 


0.3338i 


2 


2805 


+ 


0.3400i 


1 





9474 


+ 


1.0220i 


1 


6181 


+ 


1.0199i 


2 


2003 


+ 


1.0317i 


2 





5429 


+ 


1.9247i 


1 


4136 


+ 


1.7655i 


2 


0452 


+ 


1.7598i 


3 





4357 


+ 


3.1209ii 


1 


1485 


+ 


2.6098i 


1 


8314 


+ 


2.5508i 


4 





3996 


+ 


4.1798i 





8798 


+ 


3.5543i 


1 


5914 


+ 


3.4220i 


5 





3996 


+ 


4.1798« 





6275 


+ 


4.5609i 


1 


3630 


+ 


4.3662Z 


6 





3715 


+ 


5.2168« 





3095 


+ 


5.6127i 


1 


1668 


+ 


5.3579i 


7 





3502 


+ 


6.2A4U 





3160 


+ 


7.0116i 


1 


0034 


+ 


6.3729i 


8 





3335 


+ 


7.2654i 





3597 


+ 


8.0670i 





8651 


+ 


7.3974i 


9 





3202 


+ 


8.2827? 





3696 


+ 


9.1021i 





7430 


+ 


8.4249i 


10 





3093 


+ 


9.2970i 





3695 


+ 


10.129i 





6290 


+ 


9.4526ii 


11 





3002 


+ 


10.309i 





3657 


+ 


11.151i 





5140 


+ 


10.480i 


12 





2925 


+ 


11.320* 





3604 


+ 


12.169i 





3821 


+ 


11.508Z 



Table 3: QN frequencies wr/j of the tensor-gravitational perturbations for 1 = 2,1 = 3 and I = 4. 



n 




/ 




2 




/ 




3 




/ 




4 





1 


5106 


+ 


0.3575i 


2 


0079 


+ 


0.3558* 


2 


5063 


+ 


0.3550i 


1 


1 


3927 


+ 


1.1046f 


1 


9170 


+ 


1.0853* 


2 


4327 


+ 


1.0764* 


2 


1 


1938 


+ 


1.9457* 


1 


7483 


+ 


1.8697* 


2 


2915 


+ 


1.8328* 


3 





9944 


+ 


2.8990i 


1 


5380 


+ 


2.7393i 


2 


0991 


+ 


2.6477i 


4 





8460 


+ 


3.9147ii 


1 


3365 


+ 


3.6931i 


1 


8852 


+ 


3.5353i 


5 





7436 


+ 


4.9483i 


1 


1737 


+ 


4.6991* 


1 


6828 


+ 


4.4894i 


6 





6711 


+ 


5.9832* 


1 


0503 


+ 


5.7271* 


1 


5112 


+ 


5.4882i 


7 





6174 


+ 


7.0150i 





9565 


+ 


6.7613* 


1 


3730 


+ 


6.5100i 


8 





5760 


+ 


8.0431i 





8836 


+ 


7.7954* 


1 


2627 


+ 


7.5414ii 


9 





5429 


+ 


9.0678i 





8254 


+ 


8.8275* 


1 


1736 


+ 


8.5755i 


10 





5158 


+ 


10.090* 





7778 


+ 


9.8568* 


1 


1006 


+ 


9.6093* 


11 





4931 


+ 


11.109* 





7380 


+ 


10.883* 


1 


0397 


+ 


10.641* 


12 





4738 


+ 


12.126* 





7042 


+ 


11.908* 





9880 


+ 


11.671* 



frequencies, the modes are obtained up to relatively higher-order modes, associated with the mode number 
of n ~ 300. Note that the results for the scalar-gravitational perturbations are newly obtained results, 
while the results for the vector-gravitational and tensor-gravitational perturbations have been already 
calculated in j^. In the present calculations we have found no unstable mode whose frequency has a 
negative imaginary part. This is consistent with the results by Kodama and Ishibashi [22|, who proved all 
the higher-D Schwarzschild black holes to be stable against small non-radial disturbance. 

In Tables 1 through 3, the oscillation frequencies of the thirteen low-lying QNMs associated with 
/ = 2,3,4 for the scalar-gravitational, vector-gravitational, and tensor-gravitational perturbations are, 
respectively, tabulated. In the tables, the frequencies are shown in the units of the inverse of the horizon 
radius of the black hole because the horizon radius gives us one of the most natural units for the low- 
order QN frequencies. Note that another convenient frequency unit, suitable to the asymptotic behavior 
of highly damped QNMs, is the Hawking temperature of the black hole, given by Th = (n — l)/47rr/j (see, 
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Figure 1: Real parts of the non-dimensional QNM frequencies, ujvh, given as a function of the imaginary parts of 
the frequencies for the gravitational scalar-gravitational perturbations associated with 1 = 2,1 = 3, and I = 4. First 
sixty QNMs are displayed. The solid circles are used to indicate the QN frequencies. 




Figure 2: Same as Fig. 1 but for vector-gravitational perturbations. 



e.g., ini, 

In order to check our numerical code, first of all, let us compare our results with those obtained 
by Konoplya |13], who calculated fundamental QNMs in the higher-D Schwarzschild black hole with the 
sixth-order JWKB approximations. Konoplya's frequency of the QNMs associated with I = 3 for the 
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Figure 3: Same as Fig. 2 but for the tensor-gravitational perturbations. 



scalar-gravitational perturbations in the five-dimensional case is given by wr/j = 1.6044 + 0.3127i in the 
units we use, which is the only available data for the comparison. It is then found that our result is in good 
agreement with that of Konoplya, and that the relative difference between two frequencies is less than 1%. 

In the 4D Schwarzschild spacetime, there are two important properties of the QN frequencies: (i) Two 
types of gravitational perturbations, namely axial and polar perturbations, associated with the same angu- 
lar eigenvalue I possess exactly the same set of QN frequencies; (ii) There is a purely imaginary frequency. 
Interestingly, these two properties are attributed to the existence of a special relation between effective 
potentials of corresponding wave equations for scalar-gravitational and vector-gravitational perturbations. 
If we write two potentials for scalar-gravitational and vector-gravitational perturbations as V4. and V-, 
respectively, the existence of this special relation is equivalent to the existence of a function F satisfying 

V± = ±f^ + F-' + cF, (3.2) 
ar 

where c is some constant |2^. For the higher-D case, Kodama and Ishibashi found that the function F 
satisfying equation ( |3.2| ) for scalar- and vector-gravitational perturbations does not exist in the D > A 



spacetime [12|. In the five dimensional case, due to the conclusion of Kodama and Ishibashi, it is expected 
that the three decoupled types of gravitational perturbations have different sets of the QN frequencies, 
and that there is no purely imaginary frequency of the QNMs. Our numerical results are consistent with 
these predictions. On Tables 1 through 3, it is observed that the three types of gravitational perturbations 
do not have the same set of QN frequencies even when the modes are associated with the same angular 
quantum number /. We also found no purely imaginary frequency of the QNMs. 

In Figs. 1 through 3, we also display the real parts of the frequencies, Re(c(;r/i), versus the the imaginary 
parts of the frequencies, lm.{ijjrh), for the oscillation frequencies of the lower-lying QNMs associated with / = 
2, 3, 4 for the scalar-gravitational, vector-gravitational, and tensor-gravitational perturbations, respectively. 
In these figures, the first sixty QNM frequencies are exhibited. In Figs. 1-3, we can observe that the 
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real parts of the QNMs frequencies are monotonically decreasing functions of the imaginary parts of the 
frequencies except for the cases of the vector-gravitational perturbations associated with / = 3 and / = 4. 
The behavior of the QNMs of the vector-gravitational perturbations associated with / = 3 and / = 4 are 
similar to those of the gravitational perturbations in the 4D Schwarzschild black hole: the real parts of 
frequencies first decrease to a minimum value, increase to the local maximum, and then approach the 
asymptotic constant values, as the imaginary parts of the frequencies are increased. 

Finally, let us discuss the asymptotic behavior of the QNMs in the limit of highly damped modes. In 
the Schwarzschild black hole in the D = n + 2 dimensional spacetime, the asymptotic QN frequencies in 
the limit of highly damped modes are given by 

LoT-^ = log 3 + 27ri(n + ^) , (3.3) 

regardless of the spacetime dimensions and the angular quantum number I |11, |5|. On the other 



hand, higher-order corrections of the frequency ( p.3| ) depend both on the spacetime dimensions and the 
angular quantum number / |^]. For the vector-gravitational and tensor-gravitational perturbations, the 
asymptotic behavior were already investigated in detail by Cardoso, Lemos and Yoshida Q. Thus, we 
only consider the case of the scalar-gravitational perturbations, here. As shown in equation ( p. 31 ), we have 
Re(a;r;j) — > (log3)/27r as Im{ujrfi) — > oo for the 5D case. Fig. 1 shows that the real parts of the QN 
frequencies approach to the asymptotic constant value (log 3) /27r as the imaginary parts of the frequencies 
increase, and that our numerical results are consistent with the analytical prediction of the asymptotic 
behavior in the limit of large imaginary frequency. Unfortunately, we cannot however obtain the QNMs 
associated with sufficiently large mode number with our numerical code. The reason of this is unclear, but 
we may guess that this is because the recurrence relation (|2.11| ) for obtaining the QN frequencies becomes 
eight-term relation in the case of the scalar-gravitational perturbations. We therefore have failed to extract 
the detailed information of the asymptotic behavior of the QNM in the scalar-gravitational perturbations. 

4. Conclusion 

We have calculated the QNMs for the three types of gravitational perturbations of the 5D Schwarzschild 
black hole through a continued fraction method. In order to examine the QNMs, we made use of 



Schrodinger-type wave equations for determining the dynamics of the gravitational perturbations |12|. 
To apply the continued fraction method, we expanded the eigenfunctions around the black hole horizon 
in terms of a Frobenius series. It was found that scalar-gravitational perturbations obey an eight-term 
recurrence relation, and vector-gravitational and tensor-gravitational perturbations obey a four-term re- 
currence relation. For all the types of perturbations, the QNMs associated with I = 2, I = 3, and / = 4 
were calculated in the present study. Our numerical results are summarized as follows; (i) The three 
types of gravitational perturbations have diff'erent QNM frequencies; (ii) There is no purely imaginary 
frequency mode (the above two results are consistent with Kodama and Ishibashi's results that there is 
no function satisfying equation ( |3.2D ); (iii) The QNMs belonging to the three types of gravitational per- 
turbations have the same asymptotic behavior in the limit of the large imaginary frequencies, given by 
ujT^^ — > log 3 + 2'Ki{n + 1/2) as n ^ oo; This asymptotic behavior was already expected |l5[. 

For the black hole dynamics, the lower-order QNMs play an important role after the initial non-linear 
eff'ects have become insignificant. In the 4D Schwarzschild black hole, one can guess which frequency 
will be excited in the ringing phase since the axial and polar QNMs have the same frequencies. On 
the other hand, there are three different sets of QNM frequencies for the Schwarzschild black hole in 
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D > 5 dimensions. Our results show that the fundamental QNM for the scalar-gravitational perturbations 
have the smallest frequency and the longest damping-time. This means that the QNMs for the scalar- 
gravitational perturbations can live longer than other types of perturbations if all the types of perturbations 
are excited with the same initial amplitudes. Probably, in general situations things are not so simple. What 
type of perturbations will be excited will be strongly dependent on the situation. In order to expect what 
frequency will be excited in the ringing phase, one has to model the source terms of the master equations. 
So far, only the case where a test particle radially falls into a black hole has been considered |14]. In this 
case, due to the symmetry of the motion, only the scalar-gravitational perturbations are excited. The 
general case remains as a work to be solved. 
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